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ABSTRACT 

We report the first detection of CO(l— >0) emission from a submillimeter-selected galaxy, using the Green 
Bank Telescope. We identify the line in the spectrum of SMM J13 120+4242 as a broad emission feature at 
Z = 3.408, with AV FW hm = 1040 ± 190km s" 1 . If the observed CO(1->0) line profile arises from a single 
object and not several merging objects, then the CO(4^3)/CO(1^0) brightness temperature ratio of ~ 0.26 
suggests nQlj) > 3- 10 x 10 2 cm" 3 and the presence of sub-thermally excited gas. The integrated line flux 
implies a cold molecular gas mass M(H 2 ) = 1.6 x 10 M@, comparable to the dynamical mass estimate and 
4 times larger than the H2 mass predicted from the CO(4^3) line assuming a brightness temperature ratio of 
1.0. While our observations confirm that this submillimeter galaxy is massive and gas-rich, they also suggest 
that extrapolating gas masses from J upper > 3 transitions of CO leads to considerable uncertainties. We also 
report an upper limit to the mass of cold molecular gas in a second submillimeter galaxy, SMM J0943 1+4700, 
ofM(H 2 )<4x 1O 1O M . 

Subject headings: galaxies: ISM — galaxies: high redshift — galaxies: formation — galaxies: evolution — 
infrared: galaxies — radio lines: galaxies 



1. INTRODUCTION 

The resolution of the far-infrared (far-IR) background 
at long wavelengths through ground-based submillimeter- 
and millimeter-wave surveys has revealed a population of 
infrared-luminous, but optically faint, high-redshift galax- 
ies, whose spectral energy distributions (SEDs) and submil- 
limeter flux densities suggest large bolometric luminosities 
(~ 1O 13 L0), mostly emitted in infrared (IR) wave band s . Sub- 
seque nt deep X-ray observations by Alexander et al. (2003, 
2005) indicate that the bulk of this emission is not powered 
by active galactic nuclei (AGN), suggesting that high rates 
of star formation (~ 10 3 M©yr _1 ) are responsible for the ma- 
jority of the enormous luminosities of these submillimeter- 
selected galaxies (SMGs), even though they appear to already 
possess significant stellar masses (~ 1O U M ; Borys et al. 
2005). SMGs have a surface density on the sky a factor 
of 10 smaller than the optically-selected, star forming, Ly- 
man break galaxies, but are 10 times as luminous; this, to- 
gether with the f act that SMGs seem to be strongly clustered 
( Blain et al. 2004), would suggest that SMGs are an important 
population of forming galaxies, and either represent the for- 
mation epoch of some of the most massive galaxies or trace 
out unusually active regions. However, the role they play in 
the scenario of hierarchical galaxy formation and evolution 
is not yet fully understood. The large far- IR luminosity and 
significant optical obscuration of SMGs ( Smail etai]l2004l) 
are reminiscent of ultraluminous IR galaxies (ULIRGs) ob- 
served in the local universe, most of which are strong ly inter- 
acting or merging systems (Sanders & Mirabel 1996). Since 
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UL IRGs also show evidence for evolving int o elliptical galax- 
ies ( iGenzel et al ] |200"H iTacconi et al.ll2002l it has been sug- 
gested that SMGs are high-z analogs of ULIRGs, possibly 
progeni tors of the large bu lges or elliptical galaxies observed 
locally (Smail et al. 2002). However, determining the bary- 
onic and dynamical masses of SMGs, their system dynam- 
ics, the duration of the phenomenon, and their end state, and 
matching them with theoretical predictions as well as with lo- 
cal objects has proven a major challenge — the faintness of 
SMGs at optical and near-IR wavelengths precludes the de- 
tailed study needed to assess commonly-used evolutionary in- 
dicators. 

Studies of molecular gas emission provide the crucial mass 
and dynamical information needed to evaluate the evolution 
and mass of SMGs in the context of hierarchical galaxy for- 
mation and evolution. The total intensity of molecular emis- 
sion indicates the mass of gas that is available to fuel future 
star formation, and the brightness temperature ratios from ro- 
tational transitions of different quantum number J constrain 
the temperature and density of the gas, suggesting how the 
gas is being consumed. Even though millimeter-wave CO 
emission from SMGs is near currently achievable detection 



limits, a series of in terferometric surveys ( e.g. Fraveret al 
^998[fT999llNeri et al!2003tlGreve et all2005llTacconi et al 
2006) have successfully detected J > 3 CO(7 — * J- 1) tran- 
sitions in SMGs, fin ding that SMGs are massive (median 
M dyn > 7 x 10 10 M Q ; ITacconi et aTll2006l) . gas-rich (median 
M(H 2 ) = 3.0 x 1 10 M^:lGreve et alJ2005l) . and compact (me- 
dian D < 4kpc; Tacconi et al. 2006), consistent with the hy- 
potheses that SMGs are high-z counterparts of ULIRGs and 
progenitors of massive spheroids or ellipticals. All of this 
information derives from warm molecular gas, but we have 
limited knowledge of the less excited, cold molecular gas 
that makes up a substantial fraction by mass of the gas con- 
tent of present-day galaxies (e.g. the observations of M82 
of WeiB, Walter, & Scoville 2005b). With a low threshold 
for excitation, requiring n(H 2 ) ~ 10 2 cm~ 3 and AE/k ~ 5K, 
CO(1^0) emission (v res , = 115.271GHz) is the most repre- 
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sentative tracer of the metal-enriched molecular gas mass in 
galaxies be cause it is sensitive to the cold , diffuse gas that may 
dominate. Papadopo ulos & Ivisonl ( 120021) and Greve, Ivison, 
& Papadopoulos (2003) present evidence for large amounts 
of cold, low-excitation gas at high redshift, as first suggested 
by Papadonoulos et al. (2001), in the galaxy HR 10. In this 
submillimeter-bright ERO that lacks obvious AGN signatures, 
at Zco = 1-439, the CO(5^4) line luminosity underestimates 
the cold molecular gas mass by nearly an order of magnitude. 
Observations of CO(l — > 0) in SMGs are thus complementary 
to studies of / > 3 gas and necessary to determine the total 
mass of molecular gas they contain. 

Previous studies of CO(l— *0) emission from high-z galax- 
ies have been carried out with the Very Large Array (VLA) 
(e.g. Papadopoulos et al. 2000; Carilli et al. 2002b; Greve et 
al. 2003; Greve, Ivison, & Papadopoulos 2004). The Robert 
C. Byrd Green Bank Telescope (GBT) far exceeds the current 
capabilities of the VLA in instantaneous spectral bandwidth 
(800 MHz maximum per spectrometer bank, compared with 
50 MHz at the VLA) with a comparable effective collecting 
area. Since high-z CO lin es can be wider than the VLA band- 
width (Greve et al. 2004), the GBT is better suited to search 
for and measure accurately the spectral shape of CO emission 
from SMGs, at the expense of limited spatial resolution. 

In this paper, we present the results of a search for 
CO(1^0) emission from SMGs, undertaken with the GBT 
and the facility's K-band receiver (18-26 GHz). The tun- 
ing range of the GBT's K-band receiver restricts our po- 
tential SMG targets to those with redshifts greater than z — 
3.35, of which there are only four with m easured redshifts 
( Chap man etalJ 12005b iLedlow et al] 120021) . Of these four 
SMGs, we chose to target for GBT observations the two 
which also have previous high-/ CO detections (both from 
CO(4— >3) observations at the IRAM Plateau de Bure interfer- 
ometer) so we could compare the line strengths and widths 
of any detections: SMMJ0943 1+4700 at z C o = 3.346 (Neri 
et al. 2003; hereafter SMM 09431) and SMMJ13 120+4242 
at zco = 3.408 (Greve et al. 2005; hereafter SMM 13120). 
Throughout our discussion, we assume J7m = 0.3, Ha = 0.7, 
and Ho = 70 km s" 1 Mpc" 1 . 

2. OBSERVATIONS 

We searched for CO( 1 ->0) emission from SMM 0943 1 and 
SMM 13120 with the 18-26 GHz dual-feed, dual -polarization 
receiver (K-band) at the GBT on UT 2004 December 1. We 
utilized the position-switching "Nod" observing pattern, in 
which the two feeds alternate between the target position and 
a blank-sky position 3' away in azimuth such that one feed 
(beam) is always on the target source. A 2-minute integra- 
tion (scan) was completed before the receiver moved in az- 
imuth to the alternate position; thus a full Nod sequence, con- 
sisting of a pair of scans, lasted 4 minutes. Every 90 min- 
utes, the pointing and focus of the telescope were checked 
on a nearby bright radio source. Observing conditions were 
excellent through both night and day, with low atmospheric 
opacity (~ 0.02 nepers) and system temperatures in the range 
30-50 K. Total on-source integration times were 6.0 hours for 
SMM 09431 and 6.2 hours for SMM 13120. The size of the 
main beam at the GBT at 26 GHz is 28", so our observations 
are not spatially resolved. However, since 1" corresponds to 
~ 7.4 kpc at z — 3.4, it is very likely the entire extent of our 
target galaxies and any companions will be within the main 
beam (~ 2 10 kpc at z — 3.4) and our observations should be 
representative of the average cold gas properties over the ob- 



served galaxy. 

We used the GBT's Autocorrelation Spectrometer (ACS) in 
3-level sampling mode. For each of the two target sources, we 
observed a 800 MHz bandpass centered on the frequency at 
which the CO(1^0) line should fall based on the CO system 
redshift obtained from the previous CO(4^3) detection. The 
spectrometer banks were configured in low-resolution mode, 
providing 2048 channels of width 391 kHz each (~ 4.4 km s -1 
at 26 GHz) for both of the two circular polarizations of each 
feed. The noise diode for the K-band receiver was fired every 
1 .0 sec for calibration; we average the data taken with the cali- 
bration noise diode on with that taken with the noise diode off 
to improve our signal-to-noise (S/N) ratio. The noise diode 
temperature (T ca i) values for the K-band receiver provided by 
the observatory were verified using spectral observations of 
the radio-loud quasars 3C 147 (for SMM 09431) and 3C295 
(for SMM 13120) taken with the same spectrometer setup as 
the science targets. 



3. CALIBRATION 

In single-dish radio observations, the spectral baseline is the 
limiting factor in detecting very faint, broad lines like those 
expected to be produced by CO(1-^0) emission from a high- 
z galaxy. Vanden Bout, Solomon, & Maddalena (2004) dis- 
cuss the significant baseline curvature seen in wide-bandwidth 
GBT spectral data obtained over a range of elevations and 
timescales, despite the off-axis feed arm design of the tele- 
scope, which was intended to reduce reflections of signal 
along the optical path. The spectral baselines at the GBT are 
caused by a combination of small differences in the telescope 
system power between on- and off-source scans and reflec- 
tions between components in the optical system, and cannot 
be simply modeled with a low-order (n < 10) polynomial. Our 
GBT data are certainly affected by complex baseline shapes, 
since we average over ~ 6 hours of integration and 60° of el- 
evation. The calibration of our spectral data is thus nontrivial; 
we must take special care to minimize the systematic errors in 
our spectra associated with baseline ripples or curvature, and 
strive to not introduce any line-scale structure into our spec- 
trum while attempting to rem ove the baseline. 

IVand en Bout et al. (2004) outline a method for calibrating 
GBT spectral data of faint sources, in which the baseline 
shape is determined directly from the data and subtracted. 
They first produce a difference spectrum for both polariza- 
tions of each beam for every pair of scans that comprise a full 
Nod sequence according to the formula 



TdiffW) = T sys , ( 1 ) 



where ON and OFF represent an on-source scan and a blank- 
sky reference scan, respectively, observed with the same re- 
ceiver feed. The resulting spectra are averaged over time and 
feed, weighted by (7^ s )~ 2 , to produce a master average sci- 
ence spectrum T(y) for each polarization. The template con- 
taining the baseline shape of each master spectrum is con- 
structed next by first taking the scans in pairs of Nod se- 
quences (four scans) and subtracting successive ON scans and 
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FIG. 1 . — (a) Final averaged spectrum of Feed 1 for SMM 09431 data, pro- 
duced by the GBT spectral calibration method of vanden Bout et al j 120041) . 
The RMS per spectral channel is 0.37 mjy. (b) Final averaged spectrum re- 
sulting from the application of the aOFF] + /3OFF2 calibration method on 
the same data. The RMS per spectral ch annel is 0.4 1 mjy. Based on the 
CO(4^3) redshift of z = 3.346 from lNeri etaT] <2003l) . we expect to see the 
CO(1^0) line at v ohs = 26.523 GHz. 



OFF scans in the following manner: 
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where ON a and ON/, are the ON scans from the first and sec- 
ond Nod sequences, respectively, and similarly for OFF a and 
OFFb- Then, the sequence of ON-ON and OFF-OFF spec- 
tra is averaged to fo rm the master base l ine te mplate B( v) for 
each polarization. Vand en Bout et al] J2004I) point out that 
B(v) contains the same baseline structure as the science spec- 
trum T(y) (though with higher amplitude due to the longer 
time delay between differenced quantities), but without sig- 
nal from the emission line being observed. Hence, to remove 
the baseline described by the template so that the only spec- 
tral shape remaining is due to the observed source, B{v) is 
smoothed and fit to the master science spectrum for each po- 
larization by finding the best-fit coefficients (a,b,c) such that 



T(v) = aB(v) + bv+c. 



(3) 



The best-fit template is then subtracted from T(v). In Fig- 
ure 1(a) we show the final averaged spectrum after applying 
this method of baseline reconstruction and subtraction to our 
Feed 1 data on SMM 0943 1 and averaging the two polariza- 
tions of the feed together. The majority of the baseline shape 
has clearly been removed in Figure 1(a), but a slight resid- 
ual baseline remains. While this residual is less complex than 
the original baseline shape, its origin is unclear; this causes 
it to be difficult to fit and subtract, which may interfere with 
identification of faint lines from the target source. 

We preferred not to fit and subtract a baseline template from 
our final averaged spectrum at all, hoping to avoid the uncer- 
tainties associated with the process. As a result, we chose 



to use a different method of calibration, a variation of the 
traditional method of single-dish radio spectrum calibration 
described in Equation 1, which we refer to as the aOFFj + 
/3OFF2 method. In effect, we use the two OFF scans sur- 
rounding each ON scan to construct an improved, interpolated 
OFF, separately for each polarization in each feed. However, 
instead of assuming that the preceding and following OFFs 
contribute equally to the spectral shape of each ON scan, 
which might be expected if the spectral baseline was purely 
a linear function of time, we use least-squares fitting to find 
values a and (3 for each scan such that 

ON - ( a OF F x +(3 OF F 2 ) 



0, 



(4) 



aOFF 1 + (30FF 2 
where OFF\ refers to the scan prior to the ON scan and OFF 2 
refers to the scan following the ON scan (all taken through 
the same feed). Allowing a and (3 to vary by scan helps 
remove systematic shapes caused by sporadic phenomena of 
varying amplitude; we found such coefficient variation neces- 
sary because even though for ~ 70% of scans the (a, (3) values 
were (0.5 ± 0.05,0.5 ± 0.05), it was not unusual (the remain- 
ing ~ 30% of scans) for a and (3 to lie in the range 0.25-0.45 
or 0.55-0.75, where (3 = 1 -a. After the spectral fits, we apply 
the derived values of a and (3 to obtain the normalized differ- 
ence spectrum, for each scan, in units of antenna temperature: 

, s ON-(aOFF 1 + f30FF 2 ) , n s ,^ 

aOFF\ + [3OFF2 
where T sys \ is the system temperature derived for the first OFF 
scan and T sys 2 is the system temperature derived for the sec- 
ond OFF scan. The difference spectra from all scans are then 
averaged together, weighted by the root-mean-square devia- 
tion (RMS) of each scan. In Figure 1(b) we show the re- 
sult of the aOFFi + f30FF 2 method on the same spectral 
data used in Figure 1(a). The final averaged spectrum in 
Figure 1(b) lacks obvious residual ripples on the scale of a 
galaxy emission line (~ 100- lOOOkms" 1 , or - 9-90 MHz at 
v a bs = 26.5 GHz), though a long period (~ few hundred MHz), 
low amplitude ripple may be present which is not likely to 
hinder faint line detection. Note that this reasonably flat spec- 
trum was produced without subtracting a curved or rippled 
baseline template. In the case of SMM 13120, a clearly linear 
baseline appeared in the final averaged spectrum produced by 
this method, which was easily fit and removed. In general, 
we find that the aOFFi + f30FF 2 met hod produces qualita- 
tively similar results to that used by vanden Bout et al. ( 2004): 
sources detected with one calibration method are also detected 
by the other method, and the resulting spectra have compara- 
ble single-channel RMS values (the RMS values of the final 
average spectra, in the sense vanden Bout et al. method versus 
aOFFj + {30FF 2 method, are 0.37 mjy versus 0.41 mjy for 
SMM 09431 and 0.44 mjy versus 0.42mJy for SMM 13120). 
In addition, neither method can completely remove all base- 
line shape from our GBT data. Consequently, we cannot say 
that one method produces better results than the other. Even 
so, we prefer the aOFFi + f30FF 2 method because no curved 
baseline subtraction is required. However, it is important to 
note that the method is suitable only for sources with very 
weak or undetectable continu um emission at \ rest = 2.6 mm. 
SEDs of SMGs (e.g. lChapman et alJl2005l) predict that con- 
tinuum emission from SMGs at 30 GHz is undetectable in 
6 hours of integration time with the GBT. 

Before averaging the spectral data from the different polar- 
izations and feeds, we examined the periodograms (Fourier 
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power spectra) of the calibrated spectra from both polariza- 
tions and feeds of the individual scans, looking for any peaks 
that might indicate systematic spectral features not removed 
in the calibration process. Such ripples could influence the ap- 
pearance of the spectrum, and we were especially concerned 
about any ripples of width similar to the expected width of 
a CO emission line from a high-z source (~ 300km s , or ~ 
30 MHz). We found that the right polarization of feed 2 nearly 
always showed an irregularly changing pattern of power be- 
tween 0-0.03 MHz" 1 (periods > 30 MHz). This feature was 
not found in feed 1 data, and was found in approximately half 
of the left polarization data from feed 2. Thus, in an effort 
to ensure that what appears as noise on the scale of an emis- 
sion line is due to random processes, we have not included 
any feed 2 data in the averaged calibrated spectra. While this 
exclusion cuts our effective integration time in half and in- 
creases the spectral rms, it helps ensure that the quoted noise 
is uncorrected in frequency, and increases our confidence in 
the reality of any possible emission line detection by remov- 
ing ripples and structure that could obscure the line. 

To convert our object spectra to units of Janskys, we 
corrected for atmospheric attenuation and then applied the 
pa rametrization o f the GBT gain versus elevation curve found 
by Condon ( 2003 ), re-scaled for K-band by assuming an aper- 
ture efficiency at 26 GHz of 0.55 from an interpolation of 
GBT aperture efficiency measurements at nearby frequencies. 
We estimate that our flux calibration is accurate to within 
15%, due to uncertainties in the temperatures of the calibrator 
diodes, the GBT gain curve, and the calibration noise diode 
variations across the observed bandwidth, as well as the opac- 
ity estimate. 

All of the calibration of our spectral data was accomplished 
with the new library of routines written by the GBT staff in 
IDL, GBTIDL 5 . 

4. RESULTS 
4.1. SMM 13120 

SMM 13120, also known as SSA 13-332, was identified as 
a source in the Hawaii Deep Field SSA 13, with feci/im = 
6.2 mJy. The galaxy, with an optical redshift of z = 3.405, is 
not known to be lensed, and its optical sp ectrum reveals activ e 
nucleus signatures (C IV, Si IV, O III; IChapman et alJl2005l) . 
Chapman et al. (2005) estimate from the radio-submillimeter 
SED a bolometric luminosity for SMM 13120 L bot = 2.02 x 
1O 13 L and a characteristi c dust temperature of T c { = 47 K, 
though lKovacs et al.1 fe006l) find from their 350/im SH ARC- 
II observations of SMGs from the Chapman et al. ( 2005 ) sam- 
ple that the luminosities presented in Chapman et al. (2005) 
are typically overestimated by a factor of 2 and the dust 
temperatures are overestimated by 13%, which would bring 
down the values for SMM 13120 to L hol = 1.01 x 10 13 L Q 
and Td = 42 K. High-/ CO observations have confirmed 
the redshift and shown SMM 1312 to be very massive in 
molecular gas. iGreve et al.1 J2005I) detect CO(4->3) emis- 
sion from SMM 13120 with the IRAM Plateau de Bure in- 
terferometer at z = 3.408, with velocity width AVfwhm = 
530 ± 50km s" 1 and integrated line flux S co AV = 1.7 ± 
0.3 Jy kms -1 . Assuming a CO luminosity-to-H2 conversion 
factor a = 0.8M Q (Kkms~ 1 pc 2 )" 1 , this yields a molecular 
gas m ass of M(H 2 ) = (4.2 ± 0.7) x 10 10 M Q . IGreve et al.l 
J2005I) also estimate a lower limit on the dynamical mass of 

5 Information on GBTIDL is available at http://gbtidl.sourceforge.net./ 
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FIG. 2.— CO(l^O) spectrum of SMM 13120, produced by the aOFFj 
+ /3OFF2 method. A linear baseline has been subtracted from the spectrum 
post-calibration. The spectrum has been boxcar-smoothed to a velocity res- 
olution of 94 km s and re-binned so that each channel is independent. A 
Gaussian fit to the line with FWHM 1040km s is overplotted as the thick 



SMM 13120 of M dyn > 1.2 x 10 n M Q . SMM 13120 is unre- 
solved in the 6."9 x 4"8 beam at IRAM, so no information on 
the spatial structure of the molecular gas is available. 

We detect CO(1^0) emission from SMM 13120 centered 
at v bs = 26.1481GHz, with an intensity-weighted line cen- 
troid in redshift space of z = 3.408 ± 0.004. This detection 
represents the first detection of CO(l— >0) from a bona fide 
SMG and the first SMG CO(1^0) detection with spectral in- 
formation. In Figure 2 we show our calibrated and averaged 
spectrum of SMM 13120, produced by the aOFFi + /?OFF 2 
method described in §3. We have subtracted a linear base- 
line from the spectrum, fit to all the channels away from the 
line, and smoothed it from its original velocity resolution of 
4.47 km s -1 and RMS per spectral channel a = 0.42 mJy to 
94km s -1 and a = 0.16mJy. The detected line is still visible 
if the data from feed 2 are included in the final average spec- 
trum. 

The Gaussian fit overplotted in Figure 2 shows a wide 
CO(1^0) line, AV FWHM = 1040 ± 190 kms" 1 . However, the 
Gaussian fit to the line profile shown in Figure 2 appears poor, 
possibly because of the triple-peaked appearance of the line, 
although the x 2_va l ue of the Gaussian fit, \ 2 = 0.19, indi- 
cates that we cannot formally conclude whether or not the 
Gaussian function provides a good description of the line pro- 
file. To obtain a less model-dependent estimate of the line 
width, we compute also the second moment of the emission 
spectrum, finding AV = 350km s -1 . Integrating over the line, 
we find Sco^V = 0.42 ± 0.07 Jy kms -1 , where the error in- 
cludes the (estimated) 15% uncertainty in our flux calibra- 
tion. Applying the relations for Leo an d L' co presented in 
Solomon, Downes, & Radford (1992), this corresponds to a 
CO luminosity of L co = (1.0 ±0.2) x 10 7 L© oxL' co = (2.Q± 
0.3) x 10 n Kkms _1 pc 2 . Adopting the CO luminosity-to-H 2 
mass con version factor a = O.SM^ (K km s" 1 pc 2 )" 1 , which is 
found by Dow nes & Solomon! d!998t) to be appropriate for 
molecular gas exposed to the intense ultraviolet (UV) radi- 
ation fields and strong tidal forces present at the centers of 
nearby ULIRGs, we estimate the total mass of low-excitation 
H 2 in SMM 13120 to be M(H 2 ) = (1.6 ±0.3) x lO n M . Note 
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that the CO-to-H2 mass convers ion factor can v ary depend- 
ing on local gas conditions ("e.g. lYao et al.l l2003) and is thus 
highly uncertain for galaxies at high redshift. As a result, our 
absolute value of the molecular gas mass may have signifi- 
cant systematic error that is difficult to quantify. Should the 
gas in SMM 13120 be less dense and less excited than that 
of ULIRGs, the conversion factor for SMM 13120 should be 
larger and so we are underestimating the molecular gas mass; 
if the gas is more dense and more excited, the conversion fac- 
tor may be too high, and we are overestimating the gas mass. 

4.2. SMM09431 

SMM 09431 was identified by Cowie, Barger, & Kneib 
(2002) in deep submillimeter maps of the z = 0.41 galaxy 
cluster Abell 851 obtained with the SCUBA instrument on 

the James Clerk Max well Telescop e. From subsequent 

1.4 GHz observations, Ledlo wet alJ (|2002) estimated 7j = 
38 K, and suggested that two individual radio sources, iden- 
tified as H6 and H7, were associated with the SCUBA 
source. Optical follow-up observations sho wed H6 to be 
a nar row-line Seyfert 1 galaxy at z = 3.349 (Led low et al.l 
2002), while the weaker radio source H7 has a redshift 
of z = 3.3 47 from O III (T akat a et al. 2006 , in p repara- 
tion). Both lNeri et all (120031) and lTacconi et all (12006) detect 
CO(4^3) emission from H7 at z = 3.3460 ± 0.0001 with a 
line width AVfwhm = 420 ± 50km s -1 and an integrated line 
flux ScoAV =1.1 iO.Uykms -1 . H6 has not been defini- 
tively detected in CO. Assuming a CO to H2 conversion factor 
a = 0.8M© (Kkm s" 1 pcV. lNeri etaP ( 120031) infer a molec- 
ular gas mass M(H 2 ) = (2.1 ±0.2) x 1O 1O M . SMM 09431 
is lens ed by the foregroun d galaxy cluster by a modest factor 
of 1.2 ( ICowie et alJ2 002): this factor has been included in all 
stated values of mass and luminosity. 

We do not detect CO(1^0) emission from SMM 09431, 
which is apparent in the object's final calibrated spectrum 
displayed in Figure 1(b). Our observations yield a single- 
channel RMS of 0.41 mJy, which in turn corresponds to a 
statistical 3a upper limit to the integrated flux density at 
26.523 GHz of 0.056 Jy km s" 1 , assuming the CO(1^0) line 
follows the CO(4^3) line FWHM of 420 km s" 1 and that 
the line has a Gaussian profile. However, due to the un- 
certainties associated with the GBT spectral baselines and 
the flux calibration, we feel that the value of the error in 
the integrated flux, a S Av, found for SMM 09431 through ap- 
plication of the statistical error formula is probably an un- 
derestimate and thus will result in unrealistically low up- 
per limits on the CO line luminosity and molecular gas 
mass. Consequently, we estimate an alternate value of 
<tsav for SMM 09431 from the error in the integrated flux 
of SMM 13120, noting that the single-channel RMS val- 
ues of our final spectra of both objects are nearly identical 
(0.41 mJy for SMM 09431, 0.42 mJy for SMM 13120) and 
that the objects were observed in similar conditions on the 
same day at nearly the same sky frequency. We scale the 
error in S co AV of SMM 13120 {<j s av = 0.07Jykms _1 , see 
§4.1) by [AV FvraM (SMM09431)/AV fWHM (SMM 13120)] 1 / 2 
to account for the difference in CO line width between the two 
objects, and thus obtain the more conservative 3cr upper limit 
of Sco AV < 0.13Jykms _1 for SMM 09431. For the sake of 
completeness, we provide the 3cr upper limits of CO line lu- 
minosity and M(H 2 ) that result from both the statistical upper 
limit on the integrated line flux and our alternative upper limit. 
From Sco AV < 0.056 Jy km s _I , we obtain the upper limits to 
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the CO luminosity of SMM 09431 of L co < 1-1 X 10 6 L Q , or 
L' co < 2.2 x 10 10 Kkms _1 pc 2 . Assuming the CO-to-H 2 con- 
version factor a = 0.8M Q (Kkms -1 pc 2 ) -1 , the 3a upper limit 
on M(H 2 ) for SMM 09431 is then M(H 2 ) < 1.8 x 10 10 M Q . 
For the conservative estimate ScoAV < 0.13Jykms _1 , the 
corresponding upper limits are Lco < 2.5 x 1O 6 L0, L' co < 
5.1 x 10 10 Kkms _1 pc 2 , and M(H 2 ) < 4.1 x 10 10 M Q , which 
we believe to be more reasonable than the statistical upper 
limits. If we were to assume a broader CO(1^0) line, the 
upper limit will be still larger and therefore less restrictive. 

5. ANALYSIS & DISCUSSION 
5.1. CO Line Ratios and Gas Excitation 

With detections in multiple CO lines for SMM 13120, one 
of which is from CO( 1 — >0), we can use the ratio of the bright- 
ness temperatures of the different transitions to extract some 
clues about the globally-averaged conditions of its molec- 
ular gas. The task is made easier because SMM 13120 is 
not known to be strongly gravitationally lensed, so the ef- 
fects of differential lensing on the line fluxes are not likely 
to complicate our interpretation. Defining as the ratio 
of the area/velocity averaged brightness temperature of the 
CO(4^3) and CO(1^0) lines, and assuming all of the CO 
emission is from a single object, we find = 0.26 ± 0.06 for 
SMM 13120. This value of is notably lower than the ratio 
that is expected from a thermalized (i.e. relative level popu- 
lations described by principles of local thermodynamic equi- 
librium), optically thick source, t^jse = 1-0, implying that 
the 7 = 4 level of CO is not thermalized in SMM 13120. As 
a consequence, M(H 2 ) estimates for SMM 13120 calculated 
from J upper > 4 CO lines assuming the lines are fully ther- 
malized and optically thick will underestimate the mass of 
low-excitation gas in the galaxy. This is exactly what we find 
from our CO(1^0) detection of SMM 13120: our estimate 
of the molecular gas mass in §4.1, M(H 2 ) = 1.6 x 1O U M , 
is a factor of ~ 4 greater than the mass implied by the 
CO(4^3) observations if constant brightness temperature is 
assumed, M(H 2 ) = 4.2 x 10 10 M o . Moreover, the lack of a 
dominant population of thermalized high-excitation gas in 
SMM 13120 suggests that the CO luminosity-to-H 2 mass con- 
version factor we use in §4. 1 may be incorrect for this galaxy. 
iDownes & Solomonl ( 119981) derive that conversion factor for 
the central starbursts in ULIRGs, and from observ ations of the 
central IR-luminous starburst region of M82 Ce.g. lWeiBetaTl 
2005b) we expect the CO excitation in a starburst region to 
be thermalized at least to the CO 7 = 4 level. Since the 
7 = 4 level is observed to not be thermalized when aver- 
aging over SMM 1 3 1 20, the CO-to-H 2 conversion factor of 
Down es & Solomonl i 19981) is likely too low for SMM 13120, 
and the total molecular gas mass is likely to be even larger 
than the 1.6 x 1O U M we calculate in §4.1. 

Such low values of as we find in SMM 13120 are not 
unknown: in the Milky Way, ranges with position be- 
tween ~ 0.1-0.2 (Fixsen, Bennett, & Mather 1999). The 
line ratio we find for SMM 13120 is also reminiscent of 
the CO(5 -^4)/CO(l->0) rati o of 0.1 found in the ERO/SMG 
HR 10 by Greve et al. (2003), another high-z object in which 
the high-7 CO emission does not trace the low-excitation gas. 
That SMM 13120 and HR 10 are similar is somewhat surpris- 
ing, however, because HR 10 is not known to harbor an ac- 
tive nucleus, while SMM 13120 clearly does. On the other 
hand, the ratio in SMM 13120 is significantly smaller in 
comparison to those of the other high-z active galaxies with 
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CO(l^O) detections. PSS 2322+1944, APM 08279+5255, 
and 4C 60.07 are two QSOs and a radio galaxy, respectively, 
and have values of rm of 1.4, 1.5, and 0.7, indicating that 
the molecular gas they contain is nearly thermalized and that 
molecular gas masses derived from CO(4^3) and CO(l^O) 
emission will be roughly similar, which is found to be true 
dCarilli et all 12002b: Papadopoulo s et alJ l200ll IGreve etafl 
12004 . Since SMM 13120 hosts an AGN, we might expect 
its gas conditions to be more like those of the high-z active 
galaxies than galaxies lacking active nuclei; however, note 
that the active nucleus in SMM 13120 is probably signifi- 
cantly less luminous than those found in PSS 2322+1944 and 
APM 08279+5255, because, in its rest-frame UV spectrum, 
the AGN lines are weak and interstellar absorption feat ures in- 
dicati ve of a starburst are also clearly visible ( Chapm an et alJ 
2005). It is possible that the less-powerful active nucleus in 
SMM 13120 could influence the gas properties of the host 
galaxy in a different way than would a high-luminosity AGN, 
causing the global CO line ratios to differ between the SMG 
and the QSOs. It may thus be inappropriate to use the mere 
presence of an AGN as a factor in predicting global CO line 
ratios. 

From our measurement of r« in SMM 13120 we can ob- 
tain quantitative constraints on the average kinetic temper- 
ature and density of the molecular gas in this system. To 
this end, we assume the CO(1^0) and (4^3) emission re- 
sult from a single gas phas e and use a large velo city gradi- 
ent (LVG) model (based on R ichardson et alJll985l) to fit the 
observed line ratio. Data from only two CO lines generally 
do not provide enough information to place strong constraints 
on the excitation; however, since one of the lines is from the 
1 — >0 transition, we can place strong lower limits on the col- 
umn and number density of CO, and by extension, on the den- 
sity of H2. A standard LVG code, assuming a typical Galactic 
value for the [CO/H2] abundance Xco = 10~ 4 , velocity gradi- 
ent dV/dr = lkms'V 1 , and t cmb = (2.73 K)(l + z), finds 
a best-fit density for SMM 13120 of «(H 2 ) > 300 cm" 3 for 
a range of gas kinetic temperatures 7* = 42 - 54 K. This solu- 
tion is robust, such that the model finds no candidate solutions 
at other densities. However, in kinematic ally violent regions 
with intense UV radiation fields, such as we expect to find in 
ULIRGs and SMGs, the [CO/H 2 ] abundance is likely to be 
lower, so we also calculate the kinetic temperature and den- 
sity limits for Xco = 10~ 5 to cover a range of possible values. 
Leaving the assumed values of dV jdr and Tcmb unchanged, 
we obtain a slightly higher limit on the molecular gas den- 
sity of n(H 2 ) > 10 3 cm" 3 for a somewhat lower gas tempera- 
ture range of 7^ = 34-46K. Both of these sets of temperature 
and density limits point to low-density gas in SMM 13120, 
though we would guess that the solution for Xco = 10~ 5 is 
closer to reality. Also, both of the model fits result in ki- 
netic temperatures that are in good a greement with the dus t 
temperature of SMM 13120 found by Chap man et"ai] J 2005 ). 
correc ted by a factor of 13% as suggested by Kovacs et al. 
J2006T) . that is, Td = 42 K. 7* ~ Tj would be expected for an 
interstellar medium in which the gas and dust are well mixed 
and thermally coupled. Encouraged by these similar values, 
we adopt a fixed value of Tk = Td = 42 K, leave Xco I (dv jdr) 
free to vary, and once more use the LVG model to calcu- 
late the density of the molecular gas in SMM 13120, obtain- 
ing n(H 2 ) > 10 3 cm" 3 andX C0 /(dv/dr) = 1Q- 5 (kms' 1 pc~ 1 )- 1 , 
consistent with the results of the previous fits. Therefore, 
the general constraints we can place on the globally-averaged 



gas properties of SMM 13120 are n(H 2 ) > (3-10) x 10 2 cirT 3 
and Tk = 34-54K, for a realistic range of Xco I {dv jdr) = 
lO^-lO-^kms-'pc- 1 )" 1 . The mean gas density is thus likely 
to be relatively low in this system. Additionally, it is interest- 
ing to note that, if we assume the galaxy is a uniform-density 
sphere, our density limits imply a small upper limit to the ra- 
dius of SMM 13120 in the range 0.9-1.4kpc, which is con- 
sistent _wjlhjhejnedian upper limit on SMG source size that 
Tacconi et al. (2006) derive from their high-resolution images 
of high-/ CO emission from eight SMGs. 

The diffuse, warm, and low excitation gas we find 
in SMM 13120 has been observed in local ULIRGs 
( [Solomon et all 119971) and at high redshift in HR 10 
JCireve et alJ 120031) . In contrast, previous CO(l— >0) detec- 
tions of other z > 3 active galaxies have implied lower lim- 
its on the cold gas density nearly an order of magnitude 
larger. PSS 2322+1944, APM 08279+5255, and 4C 60.07 
have cold gas densities derived from r n and LVG models 
in the range n(W I > 10 3 - 10 4 cm" 3 fcarilli et all l2002bt 
Papadopou los et alJ200fl IGreve et all2004l) . 

For SMM 09431, we can use our 3rr upper limit calcula- 
tions to place only a lower limit on r^. We find > 1.2 
for ScoAV < 0.056 Jy km s" 1 , or r 43 > 0.53 for S C0 AV < 
0.13Jykms _1 . Both limits on r 43 suggest SMM 09431 has 
higher excitation, on average, than SMM 13120, though 
whether SMM 0943 1 is more similar to the high-z QSOs and 
radio galaxy with CO(1^0) detections than SMM 13120 is 
unclear. 

To examine the excitation of SMM 13120 and SMM 09431 
in another way, we plot in Figure 3 the integrated line lu- 
minosities for different CO lines (hereafter referred to as 
the CO ladder), normalized to CO(1^0), for SMM 13120 
and SMM 09431, as well as those of the other similarly- 
luminous high-z galaxies with CO(l— >0) detections. For 
SMM 09431, we use our conservative 3cr upper limit of 
Sco^V < 0.1 3Jykms~ in Figure 3. In addition, we include 
the CO data of lWeifi et al.1 (1200 5a) for SMM J16359+6612 at 
Zco = 2.517: this galaxy lacks CO(1^0) data and is an or- 
der of magnitude less luminous than the other high-z sources, 
but it has the most complete CO ladder of any SMG and thus 
provides an interesting comparison. In place of a CO(1^0) 
detection for this galaxy we use the the CO (l^O) flux im- 
plied by the LVG fit from lWeifi et alJ J2005al) . Also for com- 
parison purposes, we include in Figure 3 less-IR luminous 
(by 2-3 orders of magnitude) local galaxies commonly used 
as templates representative of quiescent galaxies (the Milky 
Way) and starbursts (the center of M82). Assuming all the 
CO emission comes from a single object, the CO ladder of 
SMM 13120 appears to rise much more slowly with increas- 
ing J than the ladders of the other high-z active galaxies, 
SMM 16359+6612, and M82's central region. The shallow 
rise of the CO ladder of SMM 13120 more resembles the CO 
ladders of the Milky Way and HR 10, and suggests excita- 
tion conditions in SMM 13120 somewhat between those of 
the inner (2.°5 < \l\ < 32.°5) and the outer (|/| > 32.°5) disk 
of the Milky Way. The shape of the CO ladder of the Milky 
Way (slow rise and turnover at lower J than the other sources) 
is caused by significant sub-thermal excitation, which is any 
molecular excitation process not governed by the principles 
of local thermodynamic equilibrium that causes the excitation 
temperature derived by comparing the population of different 
energy levels to be lower than the kinetic temperature of the 
emitting gas. While we cannot yet say that the CO ladder of 
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FIG. 3. — CO line ladder for high-z galaxies and local templates, nor- 
malized to the integrated CO(l— >0) flux. The data points for SMM 13120 
and SMM 09431 are the filled stars and the upward arrows, respectively. 
The thick dashed line represents the effect of constant brightness tem- 
perature on the CO line ratios. The CO transition data sources for 
BR 1202-0725 arelOmont eTaTHT99^) andlCarilli et ai]|2002ti) : for 4C 60.07 
Greveet al J 12 004) and Paoadoooulos et al. i2000f~for APM 08279+5255 
Downes et all 119991) an d IPaoadopoulos et alj l2TOj"fl): for PSS 2322+1944 
CarillietalJ [2 002b ) and Cox et alJ<2002J) : for HR lOIAndreani et alj EpOCl) 
and IGreve et all 120031) : for SMM 16359+6 612 IWeifi et alJ I2005J) and 
IKneibetalj f2005): for M82 M ao et alj 1200(3) and lWeiB et all I2005bl) : and 
for the Milky Wav lFixsen et all 119991) . 



SMM 1 3 1 20 turns over in a similar way as the Milky Way, the 
slow rise and low values of r^ suggest that sub-thermal exci- 
tation is important in both SMM 13120 and HR 10. The exis- 
tence of sub-thermal excitation in these objects would imply 
that observations of high-/ emission lines can provide only a 
lower limit to the molecular gas mass. 

Our findings fo r SMM 13120, along with the results of 
IGreve et all J2003I) for HR 10, highlight the need for obser- 
vations of low-/ CO emission to measure the entire molec- 
ular gas reservoir in high-z galaxies. Even though it is true 
for some high-z galaxies that / > 3 CO lines can trace low- 
excitation gas molecules as well as the highly excited gas, 
such a conclusion seems not to apply to SMM 13120 and 
HR 10. In addition, some dispersion clearly exists in the gas 
properties of galaxies hosting active nuclei at high redshift, 
since SMM 13120 has both at least a weak AGN and a signif- 
icant quantity of molecular gas not apparent in observations 
of high-/ CO transitions. 

5.2. Star Formation Timescales and Efficiency 

One of the important issues to address in understanding the 
role of SMGs in the large-scale picture of galaxy formation 
and evolution is the duration of star formation. We can use our 
new estimate of the total molecular gas mass in SMM 13120 
from CO(l— >0) to calculate the length of time the major star 
formation epoch will last from its respective star formation 
rate (SFR), assuming the SFR is constant. We define this char- 
acteristic star formation timescale as 



t SF =M(H 2 )/SFR. 



(6) 



We calculate the SF R from the fa r-IR luminosity Lfm using 
Equation 4 of lOmont etail (120011) . 

SFR = 8 MF {L FIR /\Q m LQ)M Q yx-\ (7) 

where Smf is a function of the stellar initial mass func- 
tion (IMF). We use as Lf/ R the value of Lboi found by 



Chap man etafl ( 120051) for SMM 13120, multipli ed by the cor- 
rection factor of 0.5 implied by the results of Kova cs et alJ 
J2006I) . and assume a Salpeter IMF to obtain Smf = 0.8. The 
resulting SFR and tgp we obtain, and the Lf/ R used in the 
calculation of the quantites, are listed in Table 1. Assum- 
ing all the far-IR emission is a result of star formation, we 
thus obtain a short characteristic star formation timescale 
tsF ~ 200 Myr, consistent with the brief burst needed to form 
the old spheroids and elliptical galaxies we observe at z ~ 0, 
and in agreement w ith the durat i on of SMG star formation 
bursts estimated by Sm ail et all (|2004) from modeling the 
rest-frame optical broad-band colors of SMGs ( 1 00-200 Myr). 
Even though the value of t$F we find for SMM 13120 is short, 
it is significantly longer than the star formation timescale 
calculated from the mass implied by the galaxy's CO(4-^3) 
emission and the assumption of a constant brightness temper- 
ature, which is only 50 Myr. The CO(l— >0) emission clearly 
indicates that there is more cold gas available in SMM 13120 
to fuel star formation than is suggested by observations of 
CO(4^3) and the assumption of r^ = 1. 

The continuum-to-line ratio Lfj R /L' C0 provides another 
means of investigating the evolutionary progression of 
SMM 13120, because it is commonly used to estimate the ef- 
ficiency with which molecular gas is converted to stars, in- 
dependent of the assumed value of the CO-to-H2 conversion 
factor. For SMM 1312 0, we calculate Lr-j R /L' r n using again 
the value of Lh„i fro m Chapman et al. (2005), corrected as 
suggested by K ovacs et alJ ( 120061) . for Lfi R , and we assume 
all of the CO(l— >0) emission comes from a single object. If 
we assume all the far-IR emission is a consequence of star 
formation we obtain an upper limit to the Lfi R /L' C0 ratio 
of ~ 51L Q (Kkms~ 1 pc 2 )" 1 , significantly less than the value 
(~ 190L© (Kkm s" 1 pc 2 )" 1 ) obtained with the / = 4 -> 3 line. 

In Table 1 we also list limits on tsF and Lf; R /L' C Q 
for SMM 09431 . We use the value of L FIR from 
iNeri e t al. (2003) with our CO(1^0) upper limits to esti- 
mate the continuum-to-line ratio, finding ts F < 15 Myr and 
Lfir/L' co > 680L©(Kkms~ 1 pc 2 )" 1 for the statistical up- 
per Unfit ScoAV < 0.056Jykms _1 , and t SF < 34 Myr and 
Lfir/L' co > 290Lo(Kkms~ 1 pc 2 ) -1 for the more conserva- 
tive limit Sco^V < 0.13Jykms _1 . The values found from 
our conservative limit are consistent with the results from 
CO(4^3). 

Local LIRGs and ULIRGs are known to trace a scal- 
ing relation between L' co and Lfj R in which galaxies with 
highe r IR lumino sit ies have h i gher CO luminosities (e.g. 
Youn g et all 119861) . lYao et alJ J2003I) find a median value 
ofL FIR /L' rn for LIRG s of 50 ± 30L Q (Kkm s" 1 pc 2 ) -1 , while 
Solo mon et alJ ill 9971) obtain a median value for ULIRGs of 
160 ± lSOLgtKkms^pc 2 )- 1 . The high-redshift QSOs and 
radio galaxies with CO(l— >0) detections have continuum- 
to-line ratios L F1R /L' C0 ~ 2OO-3OOL (Kkms" 1 pc 2 )" 1 from 
both the CO(1^0) and CO(4^3) lines, which is larger than 
that of SMM 13120 but comparable to SMM 09431. How- 
ever, in HR 10, L FIR /L' C0 = 130L Q (Kkms" 1 pc 2 )- 1 from 
CO(1^0), nearly an order of magnitude smaller than the 
continuum-to-line ratio calculated from its CO(5^4) line. 
Thus, the continuum-to-line ratio Lfi R /L' C0 in SMM 13120 
resembles local star-forming galaxies and HR 10, suggest- 
ing similar star formation efficiencies among them, whereas 
SMM 0943 1 appears more like the other high-z systems ob- 
served in CO( l->0). 
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It thus seems that the physical conditions in the interstel- 
lar medium of SMM 13120, while similar to another SMG 
(HR 10), are different than the high-z QSOs, 4C 60.07, and 
SMM 09431. Several possibilities could explain the differ- 
ent conditions. It is possible that the QSOs and 4C 60.07 
truly have higher star-formation efficiencies than HR 10 and 
SMM 13120. Also, the optically- or radio-luminous high-z 
galaxies and SMM 09431 could have more significant con- 
tributions to LfjR from dust heated by their active nuclei, 
causing an overestimation of the star formation efficiency 
when calculated from the continuum-to-line ratio. However, it 
should be considered that the LfiR -L' rn relation appears to be 
non-linear at high values of LpjR. iGao & SolomonH2004albl) 
suggest that, for this reason, the ratio Lfjx/L' co is not a good 
indicator of the efficiency of star formation in the most far- 
IR-bright galaxies such as high-z QSOs and SMGs. They 
argue that HCN emission is a better tracer of the dense gas 
associated with star formation and IR emission than CO due 
to its higher critical density for excitation (~ 10 5 cm~ 3 ), and 
since they find Lf/R and L' HCN are linearly correlated over 
three orders of magnitude in Lf/r, they suggest that the ratio 
Lfift/L' HCN may be a better indicator of star formation effi- 
ciency than Lfir/L' co . In this light, it may be more appropri- 
ate to use LfiR/L' HCN ratios or even Lfir/L' co ratios derived 
from J > 3 CO transitions, which trace higher-density gas 
than CO(l— »0), to compare the star formation efficiencies of 
the similarly luminous PSS 2322+1944, APM 08279+5255, 
and 4C 60.07 to SMM 13120 and HR 10. 

5.3. A Merger in SMM 13120? 

In §4.1, we found the velocity width of the Gaussian fit to 
the CO(1^0) line profile of SMM 13120 to be AV FWHM = 
1040 ± 190km s _I , which deviates sig nificantly from th e 
FWHM fitted to the CO(4->3) fine by iGreve et alJ EooBt . 
530±50kms _I . Such a discrepancy is not typically expected; 
in general, we expect the line profiles of CO transitions of 
different J will be of similar widths, from observations of 
nearby galaxies (e.g. lYao et al J 12003). However, as men- 
tioned in §4.1, our data are not of sufficient quality to be cer- 
tain that a Gaussian profile is an appropriate model for the 
line shape, so it may be misleading to compare the Gaussian 
FWHM velocity widths. To determine the reality of the factor 
of ~ 2 difference in velocity width from the fits, we instead 
directly compare the CO( 1 — > 0) spectrum of SMM 1 3 1 20 with 
the CO(4— >3) spectrum from Gre ve et all d2005 ) in Figure 4, 
where the amplitude of the 1— >0 spectrum has been multi- 
plied by a factor of 5 in both panels for ease of comparison. 
In panel (a), we compare the spectra smoothed to the same 
velocity resolution shown in IGreve et al.l (12*005). 58km s -1 , 
and it is apparent that the low S/N ratio of both spectra pre- 
clude detailed comparison. Thus, in panel (b), we further 
smooth and re-bin the line profiles of both CO transitions 
to 170kms _I , and find that the CO(4— »3) line is nearly as 
wide as the CO(l— >0) line, which suggests that the CO(4^3) 
line FWHM measu rement may have been underestimated in 
IGreve etafl (120051) . However, in panel (b) of Figure 4, we 
observe that the CO(l— >0) line extends to a lower velocity re- 
gion of the spectrum than the CO(4— >3) line, suggesting that 
there exists blueshifted (relative to the central system redshift) 
molecular gas in SMM 13120 which is less excited than the 
molecular gas at more positive velocity offsets. 

Possible explanations of this excitation variation with ve- 
locity arise from both of the dynamical scenarios discussed 
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FIG. 4. — Comparison of CO(1^0) (solid line) and CO(4^3) spectra 
(dotted line). The CO(4^3) spectrum has been provided by T. Greve. For 
clarity, the amplitude of the CO(l— >0) spectrum has been multiplied by a 
factor of 5 and the CO(4— >3) spectrum has been offset from zero by 4 mjy. 
The spectra are smoothed to a common velocity resolution of 58kms~' in 
panel (a) and 170km s in panel (b). 



bv lTacconi et alJ ll200 6) to explain the double-peaked CO line 
profiles they observe from other SMGs: a rotating disk and a 
merger of several galaxy components. If SMM 13120 were a 
single rotating disk, the blueward extension of the CO(l— >0) 
line profile could indicate a physically distinct region of cold 
gas within the disk not symmetrically distributed about the 
more highly excited region represented by the CO(4^3) 
emission, such as a single giant molecular cloud complex. On 
the other hand, in a rotating disk the rotation curve should flat- 
ten away from the center of the galaxy. In such a case, when 
the CO line profiles are observed averaged over the entire 
galaxy, the width of the CO(l— >0) profile would not increase 
over the CO(4^3) profile even if the lower excitation gas 
comes from an extended region as compared with the location 
of the CO(4^3) emission. Alternatively, if SMM 13120 is an 
interaction or a merger of several subcomponents separated 
by several hundred kms -1 in velocity space, the CO(l— >0) 
line profile could be different from the CO(4^3) profile if the 
merging components possess different excitation conditions. 

Continuum observations of SMM 13120 at different wave- 
lengths seem to support a scenario involving an interaction 
or merger of galaxy components. In Figure 5 we show previ- 
ously unpublished optical, near-infrared, and radio continuum 
images of SMM 13120, from Borys et al. (2006, in prepara- 
tion) and Fomalont et al. (2006, private communication). At 
the position of the 1 .4 GHz emission peak the HST F775 W- 
and F160W-band images show a compact, close-separation 
(~ 0"2) double source. Plus, the radio continuum emission 
appears to be extended on the scale of ~ 3" (at PA. ~ 42°), 
significant at the 4er level. In fact, the northern source ~ 2. "5 
away in the HST images appears to be associated with this ex- 
tended radio emission. The radio continuum indicates signifi- 
cant activity (from black hole accretion or star formation) out- 
side the radio/optical core of SMM 13120, and is suggestive 
of a galaxy interaction or merger. However, it is not necessar- 
ily associated with the blueward-extending CO(l— ►()), since 
the activity causing the extended radio emission would likely 
cause increased CO excitation. Obtaining high-accuracy red- 
shifts of the optical sources around the optical/radio core 
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FIG. 5.— 8" X 8" images centered on the radio position of SMM 13120 (0,2000 = 13 ,! 12'"01? 172, 5, 2 ooo = 42°42'08;'39). Left Panel: 8" X 8" HST-ACS 
F775W-band (i') image (Borys et al. 2006, in preparation). Center Panel: HST-NICMOS F160W-band image (Borys et al. 2006, in preparation). Right Panel: 
VLA 1.4 GHz contours plotted over a Subaru-SUPRIMECam R-band image (Fomalont et al. 2006, private communication). The VLA data were obtained in 
A-array, with synthesized beam size 1"5, and the contour levels in the image are 15 and 30^tJy. 



and/or interferometric imaging at high spatial resolution of 
the CO emission lines can verify the merger/interaction sce- 
nario: if the optical/radio sources are found to have redshifts 
corresponding to the center and outer regions of the CO( 1 — »0) 
and CO(4— >3) lines, then the sources may indeed be interact- 
ing galaxies. 

More data are clearly needed to determine the structure 
and energetics of the major gas component(s) in SMM 13120, 
and SMGs as a class. In particular, spatially- and spectrally- 
resolved images and spectra of CO(1^0) and high-/ emission 
will be the most revealing, along with HCN observations, and 
new instruments and facilities planned for the future will help 
provide these. In the near future, the new Ka-band correla- 
tion receiver at GBT will cover the frequency range to which 
CO(l— >0) is redshifted for galaxies in the range 2 < z < 3 
(26-40 GHz) and thus will aid the detection of CO(l— >0) in 
many more SMGs in a modest integration time (~ 6 hours or 
less). The GBT will also be important in obtaining HCN data 
for SMM 13120 and SMGs in general, in a few years assisted 
by the EVLA 6 . Looking further in the future, ALMA 7 will 
be instrumental in understanding the nature of SMM 13120 
and other SMGs, since it will provide spatial resolution com- 
parable to HST at submillimeter wavelengths and allow us to 
determine if and where there are multiple components con- 
tributing to the far-IR continuum emission. ALMA will also 
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tions can be explored. 
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TABLE 1 

Physical Parameters of Observed SMGs 



Source 


zco 




Lfir 
(1O 13 L ) 


SFR b 

(Moyr 1 ) 




CO Trans 


M(H 2 ) 
(1O 1O M ) 


(Myr) 


Lfir/L^o 
(Lq (Kkms-'pc 2 )- 1 ) 


Reference 


SMM 13120 


3.408 


1.0 


1.0 


810 


0.26 ±0.06 


1^0 


16±3 


200 


51 


d),(2),(3) 














4^3 


4.2 ±0.7 


50 


190 




SMM09431 d 


3.346 


1.2 


1.5 


1200 


> 0.53 


l-»0 


<4.1 


<34 


>290 


(D,(4) 














4^3 


2.1 ±0.2 


18 


560 





References. — (1) Greve el al. 2005; (2) Chapman el al. (2005); (3) Kovacs el al. 2006; (4) Neri el al. 2003 



Gravitational lensing magnification factor applied to Lpjg and Lqq, and thus included in all quantities calculated from those luminosities. 
Derived from Lp jjf according to Equation 7 in the text. 
Calculated using Equation 6 in the text. 

"Conservative" CO(l — »0) limits from §4.2 used in all calculations. 



